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Abstract 



In the framework of potential models for heavy quarkonium, we compute the 
mass spectrum of the bottom-charmed B c meson system and spin-dependent 
splittings from the Schrodinger equation using the shifted-large-N expansion 
technique. The masses of the lightest vector B+ and pseudoscalar B c states 
as well as the higher states below the threshold are estimated. Our predicted 
result for the ground state energy is 6253^g 5 MeV and are generally in ex- 
act agreement with earlier calculations. The parameters of each potential 
are adjusted to obtain best agreement with the experimental spin-averaged 
data (SAD). Our results are compared with the observed data and with the 
numerical results obtained by other numerical methods. 



I. INTRODUCTION 

Recently, theoretical interest has risen in the study of the B c meson, the heavy cb quarkonium 
system with open charm and bottom quarks composed of two nonrelativistic heavy quarks. 
The spectrum and properties of the cb systems have been calculated various times in the 
past in the framework of heavy quarkonium theory [1]. Moreover, the recent discovery of the 
lightest vector meson [2] has inspired new theoretical interest in the subject [3-6]. For the 
mass of the lightest vector meson, the CDF Collaboration quotes M Br = 6.40 ± 0.39 ± 0.13 
GeV. This state should be one of a number of states lying below the threshold for emission 
of B and D mesons. Such states are very stable in comparison with their counterparts in 
cc and bb systems. A particularly interesting quantity should be the hyperfine splitting that 
as for cc case seems to be sensitive to relativistic and subleading corrections in a s . For the 
above reasons it seems worthwhile to give a detailed account of the Schrodinger energies 
for cc, bb and cb meson systems below the continuum threshold. Because of the success 
of the nonrelativistic potential model and the flavour independence of the q{q 2 potential, 
we choose a set of phenomenological and a QCD-motivated potentials. We insist upon 
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strict flavor-independence of its parameters. We also use a potential model that includes 
running coupling constant effects in both the spherically symmetric potential and the spin- 
dependent potentials to give a simultaneous account of the properties of the cc, bb and cb 
systems. Since one would expect the average values of the momentum transfer in the various 
quark- ant iquark states to be different, some variation in the values of the strong coupling 
constant and the normalization scale in the spin-dependent potential should be expected. 

We have made this study mostly a fully treatment for the potentials used in the literature. 
In order to minimize the role of flavor-dependence, we use the same values for the coupling 
constant and the renormalization scale for each of the levels in a given system and require 
these values be consistent with a universal QCD scale. 

In 1991 Kwong and Rosner [7] predicted the masses of the lowest vector (triplet) and 
pseudoscalar (singlet) states of the B c systems using an empirical mass formula and a loga- 
rithmic potential. Eichten and Quigg [1] gave a more comprehensive account of the energies 
and decays of the B c system that was based on the QCD-motivated potential of Buchmuller 
and Tye [8] . Gershtein et al. [9] also published a detailed account of the energies and de- 
cays of the B c system and used a QCD sum-rule calculations. Baldicchi and Prosperi [6] 
have computed the cb spectrum based on an effective mass operator with full relativistic 
kinematics. They have also fitted the entire light-heavy quarkonium spectrum. Fulcher [4] 
extended the treatment of the spin-dependent potentials to the full radiative one-loop level 
and thus included effects of the running coupling constant in these potentials. He also used 
the renormalization scheme developed by Gupta and Radford [10] . 

One of the important goals of the present work is to extend the shifted large-N expansion 
technique (SLNET) developed for the Schrodinger equation [11, 12] to reproduce the cb 
spectroscopy using a class of three static together with Martin and Logarithmic potentials 
which was already utilized for producing the spin-averaged data (SAD) self-conjugate (qq) 
mesons and the (SAD) non-self conjugate (qiq 2 ) mesons [13] . We also extend our work using 
an improved QCD-motivated potential proposed by Buchmuller and Tye [8] . 

The contents of this article is as follows: in section II, we present the solution of the 
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Schrodinger equation using the SLNET for the the non-self conjugate q{q 2 mass spectrum. 
In section III we present all the potentials used in this work. A brief discussion and conclusion 
appear in section IV. 



II. WAVE EQUATION 

In this section we shall consider the N-dimensional space Schrodinger equation for any 
spherically symmetric potential V(r). If ^(r) denotes the Schrodinger's wave function, a 
separation of variables ^(r) = Y e>m (9, 0)-u(r)/r( 7V_1 ) //2 gives the following radial equation (in 
units fr=l)[ll,12] 

f 1 d 2 [k - (1 - a)][k - (3 - a)} ~\ r { \ \ { \ r, { \ 

where [i = {m qi m q2 ) / (m qi +m q2 ) is the reduced mass for the two interacting particles. Here 
E n ^ denotes the Schrodinger binding energy, and k = N + 2£ — a, with a representing a 
proper shift to be calculated later on and £ is the angular quantum number. We follow the 
shifted 1/N or 1/k expansion method [12 ] by defining 
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and also the energy eigenvalue expansion [12,13] 

E n , t =E + Ejk + E 2 /t + E 3 /t + O (l/fc 4 ) , (3) 

—1/2 

where x = k (r/r — 1) with r is an arbitrary point where the taylor expansions is being 
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performed about and Q is a scale to be set equal to k at the end of calculations. Inserting 
equations (2) and (3) into equation (1) yields 

1 d 2 1 (k (2-o) (l-a)(3-q) 



2fidx 2 2fi V 4 2 4 k 



2x 3x 2 

X 1-T- + 



fcV2 k 



<Pn r {x) = £n r <Pn r {x), 



(4) 



Q V " a:" 2 2/,- 

where the final analytic expression for the 1/k expansion of the energy eigenvalues appro- 
priate to the Schrodinger particle is [11] 



r\k 

~Q 



_, E-\ Eo E% „ 



(5) 



fc A; 2 k 3 

where n r is the radial quantum number. Here, we formulate the SLNET (expansion as 1/k) 
for the nonrelativistic motion of spinless particle bound in spherically symmetric potential 
V(r). The resulting eigenvalue of the N-dimensional Schrodinger equation is written as [12] 



in r — k 



1 rgK(rp) 
8/i Q 
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where a^and a^are the expressions given by Imbo et al [11]. Comparing equation (5) with 
(6) yields 

Q 



Eo = V(r ) 



8/ir 2 
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E 3 = —a 
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Here the quantity r is chosen so as to minimize the leading term, E [12,13]. That is, 

= and : n > 0. (11) 



dr 

Therefore, r satisfies the relation 



drl 



g = 4/ir 3 \/'(r ), 



(12) 



and to solve for the shifting parameter a, the next contribution to the energy eigenvalue E\ is 
chosen to vanish [11] so that the second- and third-order corrections are very small compared 
with the leading term contribution. The energy states are calculated by considering the 
leading term E and the second-order and third-order corrections, it implies that 



a = 2 — 2(2n r + 1)[mv, 



(13) 



with 



uo = 



4/x 



3 + 



roV'{r ) 



nl/2 
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V"(r ) 

Once r is being determined, with the choice k = \fQ which rescales the potential, we derive 
an analytic expression that satisfies r in equations (12)-(14). We have defined the analytic 
expression for the 1/k expansion of the energy eigenvalues (3) and determined the shifting 
parameter requiring E\ — [12, 13]. For the Coulomb potential, considered as a testing 
case, the results are found to be strongly convergent and highly accurate. The calculations 
of the energy eigenvalues were carried out up to the second order correction. Therefore, the 
Schrodinger binding energy to the third order is 
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Once the problem is collapsed to its actual dimension JV = 3, it simply rests the task of 
relating the coefficients of our equation to the one-dimensional anharmonic oscillator in order 
to read the energy spectrum. For the N = 3 physical space, the equation (1) restores its 
three-dimensional form, one obtains 



1 + 2£ + (2n r + 1) 



3 + 



r V'(r ) 



1/2 
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(16) 
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Once r is being determined through equation (16), the Schrodinger binding energy of the 
q{q 2 system in equation (15) becomes relatively simple and straightforward. We finally 
obtain the total Schrodinger mass binding energy for spinless particles as 



M{qiq 2 ) = m qi + m g2 + 2E nA 



(17) 



where n = n r + 1 is the principal quantum number. As stated before in [11-13], for a fixed 
n the computed energies become more accurate as £ increases. This is expected since the 
expansion parameter 1/N or 1/k becomes smaller as £ becomes larger since the parameter 
k is proportional to n and appears in the denominator in higher-order correction. 



III. SOME POTENTIAL MODELS 



A. Static potentials 

The cb system that we investigate in the present work is often considered as nonrela- 
tivistic system, and consequently our treatment is based upon Schrodinger equation with a 
Hamiltonian 

H = ~ + V(r) + V SD , (18) 

where Vsd is the spin-dependent term having the simple form 

V SD — V SS = -^^5 3 (r) Sl .s 2 . (19) 
9m qi m q2 

The spin dependent potential is simply a spin-spin part and this would enable us to make 
some preliminary calculations of the energies of the lowest two S-states of the cb system. 
The potential parameters in this section are all strictly flavor-independent. The potential 
parameters are fitted to the low-lying energy levels of cc and bb systems. The strong cou- 
pling constant a s is fitted to the observed charmonium hyperfine splitting of 117 MeV. The 
numerical value of a s is dependent on the potential form and found to be compatible to 
the other measurements of reference [1,4,6-7,9]. The hyperfine observed splitting in charmo- 
nium fixes a s for each potential. The perturbative part of such a quantity was evaluated at 
the lowest order in a s . Baldicchi and Prosperi [6] used the standard running QCD coupling 
expression 
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a,(Q) = — j-^-. (20) 

(11 In (£) 

with rif = 4 and A = 0.2 GeV cut at a maximum value a s (0) = 0.35, to give the right 
J/ip — rjc splitting and to treat properly the infrared region. Details on their numerical work 
are given in Ref. [6]. Whereas Brambilla and Vairo [3] took in their perturbative analysis 
0.26 < a s (n = 2GeV) < 0.30. 

The central potentials in (18) include a class of a static potentials of the general form 

V(r) = -ar-P + br p + c, < (5 < 1, (21) 

previously proposed by Lichtenberg [14] where the parameters a and b are positive while c 
may be of either sign. These static quarkonium potentials are monotone nondecreasing, and 
concave functions which satisfy the condition 

V'(r) > and V"(r) < 0. (22) 

This comprise the following potentials: 



1. Cornell potential 

The QCD-motivated Coulomb-plus-linear potential (Cornell potential) [15] 

V c (r) = -- + br + c, (23) 
r 

with the adjustable parameters as 

[a, b, c] = [0.52, 0.1756 GeV 2 , -0.8578 GeV]. (24) 

Such a potential has the Coulomb-like behavior at low distance is supposed to represent 
the short range gluon exchange and the term confining the quarks rises linearly at large 
distances is supposed to represent the confinement string tension. The main drawback of 
this potential is that the cc and bb states lie in an intermediate region of quark separation 
where neither limiting forms of (22) should be valid. 
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2. Song-Lin potential 

This phenomenological potential was proposed by Song and Lin [16] having the form 

V SL {r) = -ar- 1/2 + br 1/2 + c, (25) 
with the adjustable parameters 

[a, 6, c] = [0.923 GeV 1/2 , 0.511 GeV 3/2 , c = -0.798 GeV] . (26) 
The characteristic feature of this potential may be traced in Ref. [16]. 

3. Turin potential 

Lichtenberg and his co-workers [17] suggested such a potential which is an intermediate 
between the Cornell and Song-Lin potentials. This potential [17] has the form 

V T (r) = -ar- 3 / 4 + br 3 ' A + c, (27) 

with its adjustable parameters 

[a, b, c] = [0.620 GeV 1 ^, 0.304 GeV 7/ \ -0.823 GeV] . (28) 

4- Martin potential 
The phenomenological power-law potential of the form 

V M {r) =b{r xlGeV) - 1 + c, (29) 
is labeled as Martin's potential [18] with the parameters values 

[6, c] = [6.898 GeV 11 , -8.093 GeV] . (30) 
(potential units are also in GeV). 



5. Logarithmic potential 



A Martin's power-law potential reduces [19] into a form 



V L (r) =Mn(r x IGeV) + c, 



(31) 



with 



[b, c] = [0.733 Ge^, -0.6631 GeV] . 



(32) 



The potential forms in (29), and (31) were used by Eichten and Quigg and all of these 
potential forms were also used for ip and T data probing 0.1 fm < r < 1 fm region. 



The Coulomb-plus-linear potential has been successful in describing the spectrum of 
existing quarkonium systems. Moreover, the same is true for logarithmic potential as well as 
for power law potentials [8]. These potentials all describe well the shape of the intermediate 
region. Neither the long-distance nor the short-distance nature of the interquark potential 
is probed by upsilon and charmonium systems. Other systems depend on the short- and 
long- distance parts of the potential as in toponium system. The potentials for some systems 
should account properly for the running a s at the different distances, the a s and the string 
constant that describe the long-distance behavior can vary independently. Furthermore, the 
Aqcd used to evaluate a s is related to a specific renormalization scheme so that comparison 
with other calculations is possible. This leads to expanding a s to at least 2-loop order and 
use of nontrivial methods in the interpolation between long-and short-distance behaviors. 
The interquark potential at short distances has been calculated to 2-loop calculations in the 



B. QCD-motivated potentials 



1. Igi-Ono potential 
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modified minimal-subtraction (MS) scheme [20]. Together with the 2-loop expression for 
a s , one has [8,21] 



3r 



l a a (fi) 



H — - — (b lrifj,r + A) 
Ztc 



(33) 



and 



a s (n) 



4tt 



b ln{^/K?—) 



b-y In /n(/i 2 /A- 



MS' 



(34) 



with 



[6d,Ml 



2 38 , 31 5 

11 - -n /? 102 - —n f , b QlE + — - -n f 



(35) 



Here n/ is the number of flavors with mass below /x and 7^ = 0.5772 is the Euler's number. 
The renormalization scale \i is usually chosen to be 1/r to obtain a simple form for V(r). 
The singularity in a s (l/r) at r=l/A I ^ = 5 GeV' 1 m 1 fm for = 200 MeV. This 
singularity can be removed by the substitution 



In 



1 



r 2 A- 

MS 



f(r) = In 



r 2 A- 

MS 



(36) 



The constant b is an adjustable parameter of the potential and will not affect the perturbative 
part of the potential. Hence, setting rif = 4 in equation (35), the one-gloun exchange part 
of the interquark potential simply takes the form 



V, 



(n /= 4) 



OGE 



r) 



16tt 1 
25 rf(r) 



53 1 



l 462 In f{r) | 2 lE + ff 



(37) 



625 /(r) /(r) 

where the function /(r) can be read off from (36). Furthermore, the long distance interquark 
potential grows linearly leading to confinement as 



V L (r) = ar. 



Igi and Ono [21] proposed a potential whose general form 



(38) 



\/W= 4 )(r) = V^- 4) +ar + dre~ 9r , 



(39) 
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so as to interpolate smoothly between the two parts. They added phenomenologically a 
term dre~ gr to the potential so that to adjust the intermediate range behavior by which the 
range of is extended keeping linearly rising confining potential. Hence, the runs 
from 100 to 500 MeV keeping a good fit to the cc and bb spectra. Thereby, the potential 
with b = 20 is labeled as type I, the one with b = 5 is labeled as type II, and that one with 
d = and b = 19 is labeled as type III. 



2. Improved Chen-Kuang potential 



Chen and Kuang [22] proposed two improved potential models so that the parameters 
therein all vary explicitly with h-jj-g so that these parameters can only be given numerically 
for several values of A-j^s- Such potentials have the natural QCD interpretation and explicit 
Ajfg dependence both for giving clear link between QCD and experiment and for convenience 
in practical calculation for a given value of Aj^. It has the general form 

16tt 1 



V^ =i \r) = kr- 



25 rf{r) 



1 _ 462 lnf{r) + 2 lE + ff 



(40) 



625 f(r) f(r) 

where the string tension is related to Regge slope by k = The function f(r) in (40) can 
be read off from 



f(r) = In 



A 



+ 4.62 - A(r) 



MS 1 



and 



A{r) 

with parameters values 
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MS 



4Ai 



MS 
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MS 



A 



MS' 



(41) 



(42) 



[fc,a',Ay = [0.1491 GeV 2 , 1.067 GeV' 2 , 180 MeV] 



(43) 



The details of this potential can be tracesd in Ref. [22]. 
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IV. DISCUSSION AND CONCLUSION 



We solve the Schrodinger equation for various potentials to determine the position of the 
IS spin-averaged data (SAD) for cc, bb, and then for cb systems. The hyperfine splitting 
of the ground state is given by (19). The hyperfine splitting observed in the charmonium 
family [23] 

M(J/iP) - M(r] c ) = lllMeV, (44) 

fixes the strong coupling constant a s for each potential. For simplicity and for the sake 
of comparison we neglect the variation of a s with momentum in (20) to have a common 
spectra for all states and scale the splitting of cb and bb from the charmonium value in (44). 
The effective constant value, fixed by the described way, has been applied to the description 
of not only the cc system, but also the bb, and possibly cb quarkonia. The consideration 
of the variation of the effective Coulomb interaction constant becomes especially essential 
for the T particle, for which a s (T) ^ a s (ip). Thus, in calculating the splittings of the cb 
spectra, we have to take into account the a s {ji) dependence on the reduced mass of the heavy 
quarkonium instead of a s (Q) for the reasons stated in Ref. [3]. That is, the QCD coupling 
constant a s in (19)-(20) is defined in the Gupta-Radford (GR) renormalization scheme [10] 

Q 

a s = 7 v ' ( 45 ) 

(33 - 2n f ) In (^) 

in which Aqr is related to Ajjg by 

"49 - lOnj/3 



A Gj r = Ajgrrexp 



(46) 



2(33-2n/)_ 

The treatment of our model with momentum dependence form (20) would increase the 
accuracy, it probably give very close results and might reproduce the experimental values 
equally well within the errors. 

Table I reports our prediction for the Schrodinger mass spectrum of the four lowest cb 
S— states together with the first three P— and D— states for various potentials. If we use 
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the a s determined for the cc and bb systems by other authors [1,4,6-7,9], we predict their 
energy masses to within a few MeV to the calculated SAD [13,17,24]. Since the model is 
spin independent and as the energies of the singlet states of quarkonium families have not 
been measured, a theoretical estimates of these unknown levels introduces uncertainty into 
the calculated SAD. Its worthwhile to note that SAD is defined as the average mass of the 
(s = l,£=l) states in the form SAD(nPj) = § [5M(n 3 P 2 + 3M(n 3 Pi + M(n 3 P )] and for 
(s = 1J = 0) states by SAD(nSj) = \ [3M(n 3 Si) + M(n%)] in which the SAD S level 
gives the weight of only 1/4 to the unknown singlet level and 3/4 to the known triplet level. 

Instead of showing the calculated masses in GeV as in Ref. [6], it is useful to report 
the spectrum for heavy meson masses in MeV obtained by our formalism in numerical 
form to four significant figures. Obviously, this trend provides a measure of the accuracy in 
reproducing the experimental and the calculated SAD as one can expect in this match to 
a nearest few MeV . It has been shown the possibility of producing cb mesons in e + e~ and 
hadron-hadron colliders [25,26], so that the study of the cb mesons is not merely academic. 
It is clear that the differences between the spectra predicted by different potentials are not 
large. 

Moreover, the calculated fine and hyperfine splitting values of the vector and pseudoscalar 
masses of the cb system for the two-lowest S— states are also presented in Table I. Predictions 
for the cb ground-state masses depend little on the potential. The B c and B* masses and 
splittings lie within the ranges quoted by Kwong and Rosner [7] in their survey of techniques 
for estimating the masses of the cb ground state. Thus our results in Table I are very similar 
to those presented in Table I of Eichten and Quigg [1] and also to that presented by Fulcher 
[4]. The results obtained with the Song-Lin and the Martin potential in all cases fall between 
the extremes defined by the other potentials. Here, we report the range of the coupling 
constant we take in our analysis 0.220 < a s < 0.313 for all potentials and 0.264 < a s < 0.313 
for the class of static potentials in our study as shown in Table I. Overmore, our predictions 
to the cb masses of the lowest S-wave (singlet and triplet) and also their splittings appear 
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together with those estimated using various methods by many authors are also shown in 
Table II. Larger discrepancies among the various methods occur for the the ground and 
excited states [6]. As noted in Table II, Bambilla and Vairo [3] calculated the maximum 
final result of (M B *) peH = 6326^ 9 MeV, the upper limit corresponds to the choice of 
parameters A^g 3 = 350 MeV and /i = 1.2 GeV, while the lower limit to AgJ 3 = 250 MeV 
and \i = 2.0 GeV as the best approximation to their perturbative calculation. Its worthwile 
to note from Table II, that the SLNET estimations fall in the range demonstrated by other 
authors. 

The fitted parameters of the potential of type I and II which were proposed by Igi and 
Ono are listed in Table III. Our predictions for the cb mass spectra are reported in Table IV 
for the type I (with b = 20) together with the type II (with 6 = 5). Moreover, the singlet 
and triplet masses together with the hyperfine splittings predicted for the two types of this 
potential are also reported in Table V. We, hereby, tested acceptable parameters for A^ 
from 100 to 500 MeV for the type I and type II Igi-Ono potentials to produce the cb masses 
and their splittings. Small discrepancies between our prediction and SAD experiment [24] 
can be seen for higher states and such discrepancies are probably seen for any potential 
model and it might be related to the threshold effects or quark-gloun mixings. The fitted 
set of parameters for the Igi-Ono potential (type III ) presented by Ref. [21] are also tested 
in our method with b = 19 and (A^ = 300 MeV and also 390 MeV) and then b = 16.3 and 
A-jig = 300 MeV which seems to be more convenient than A^ = 500 MeV used by other 
authors [8]. Results of this study are also presented in Table VI. It is clear that the overall 
study seems likely to be good and the reproduced masses of states are also reasonable. Once 
the experimental masses of the cb spectra becomes available one may be able to sharpen 
analysis by adjusting potential parameters and then choosing the mostly convenient value 
of Aj^. We see that the quark masses m c and are sensitive to the variation of Ajfg, their 
values increase and explained in the following way. The absolute value of the short-range 
potential decreases with increasing A-jjg. In order to reproduce experimental masses we need 
larger quark masses for larger values of Aj^. This situation becomes completely different if 
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we are allowed to add an additional constant to (39) and reduce the effect of the exponential 
term therein for the higher states. However, we have not also attempted such possibility in 
this paper for the sake of comparison with the other authors. 

We have found the q{q 2 potentials which reproduce the experimental masses of the cb 
states for various values of A.jfg. Using this model, we see that the experimental cb splittings 
can be reproduced for ~ 200 - 400 MeV and b = 20 while = 100 and 500 MeV 
are clearly ruled out. Taking the preferred value of A^ = 400 MeV, we can predict the 
splittings in exact agreement to several MeV with the other formalisms (c.f., Table I of 
Ref. [6]). We have selected some preferred values for Aj^g and b as they provide the best 
close fine and hyperfine splittings to the other works as shown in Table VI. 

The predicted cb spectra obtained from the Chen-Kuang potential are also listed in Table 
VI. We see that for states below the threshold, the deviations of the predicted spectra from 
the experimental SAD are within several MeV. We find that m c and mj, are insensitive to 
the variation of Aj^g for this potential. This is consistent with the conventional idea that, 
for heavy quarks, the constituent quark mass is close to the current quark mass which is 
A^r independent. The Ajfg dependence of the Chen-Kuang potential is given in (41)- (42) 
and the predicted spectra is obtained for various values of A^ = 100, 180 and 375 MeV. 
Numerical calculations show that this potential is insensitive to Ajjg in the range from 100 
to 375 MeV, and as A^ runs from 375 to 500 MeV,the potential becomes more sensetive. 
The obtained n 1 ^ and n 3 Si splittings for the cb system in the Chen-Kuang potential are 
also listed in Table VI. For a certain range of Ajjg the agreement is good. Moreover, the 
fine and hyperfine splittings of the S-states in the cb system predicted by Igi-Ono and Chen- 
kuang potentials are listed in Table VII for some proper parameters. They are considerably 
smaller than the corresponding values A ls (cF) = 76 MeV, and A 2 s(cF) = 42 MeV predicted 
by the quadratic formalism. Moreover, Chen-Kuang [22] predicted A is (cb) = 49.9 MeV, and 
A 2 s(c&) = 29 A MeV for their potential with A^ = 200 MeV in which the last splitting is 
almost constant as A^r increases. Our predictions for Ais(cF) = 68 MeV, and A 2 s(c6) = 35 
MeV for the Chen-Kuang potential with A^ runs from 100 into 375 MeV. We also find 
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A is (cb) = 67 MeV, and A 2S {cb) = 33 MeV for the Igi-Ono potential with = 300 
MeV and b = 16.3. The present model has the following features: (1) The present potential 

predicts smaller Aig and A 2 g than the other potentials do for cb system and the present Ais 
and A 2 s do not depend on more sensitively (2) The experimental cb spilitting can be 

repoduced for the preferred range of A^r runs from 100 into 375 MeV . Table VII reports 
our results using SLNET compared to other formalisms. 

In this paper, we have developed the SLNET in the treatment of the cb system using 
group of static and QCD-motivated potentials. For such potentials the method looks quite 
attractive as it yields highly accurate results. The convergence of this method seems to be 
very fast as the higher corrections to energy have lower contribution. It is interesting to 
note that the scope of this method can be extended to more realistic potentials. 

In this respect, in reproducing the SAD, we used the same fitted parameters of the other 
authors for the sake of comparison and also for the sake of testing the accuracy of our 
approach. Here, we would expect much better agreement to experimental data in case of 
fitting our own parameters properly. Finally, we comment that the fitted parameters of any 
potential are model-dependent in any study. 
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TABLE I. The cb masses and hyperfine splittings (A nS )* 
calculated for some potentials (all in MeV). 



States [6,24] Cornell Song-Lin Turin Martin Logarithmic 



Oi s = 




0.313 


0.264 


0.286 


0.251 


0.220 


m c (GeV) = 




1.840 


1.820 


1.790 


1.800 


1.500 


m b (GeV) = 




5.232 


5.199 


5.171 


5.174 


4.905 


M(cb) 














IS 


6315 


6315 


6306 


6307 


6301 


6317 




6334 


6335 


6325 


6326 


6319 


6334 


1% 


6258 


6253 


6248 


6249 


6247 


6266 


A 1S * 


77 


82 


76 


77 


72 


68 


2S 


6873 


6888 


6875 


6880 


6892 


6903 


2 3 S 1 


6883 


6897 


6884 


6889 


6902 


6911 


2'So 


6841 


6860 


6850 


6852 


6865 


6879 


A 2S 


42 


37 


34 


36 


37 


31 


3S 


7246 


7271 


7209 


7246 


7236 


7225 


A O 

4o 




7587 


7455 


7535 


T A O O 

7483 


7448 


IP 


6772 


6743 


6733 


6731 


6730 


6754 


2P 


7154 


7138 


7104 


7123 


7125 


7127 


3P 




7464 


7371 


7428 


7398 


7375 


ID 


7043 


7003 


6998 


6998 


7011 


7027 


2D 


7367 


7340 


7284 


7320 


7311 


7301 


3D 




7636 


7510 


7588 


7536 


7502 


*A nS = M(n 


3 Si) - 


- M(n l 


So). 
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TABLE II. The predicted cb masses of the lowest S-wave and its splitting 

compared with the other authors (all in MeV). 

Work Mb^SqY M b «(1 3 S 1 ) A 1S 



Eichten et al. [1] 6258 ± 20 

Colangelo and Fazio [3] 6280 
Baker et al. [27] 6287 
Roncaglia et al. [27] 
Godfrey et al. [1] 6270 
Bagan et al. [1,27] 6255±20 
Bambilla et al. [3] 
Baldicchi et al.[6] 
SLNETt 
SLNET* 

SLNET 6310 T 



6350 
6372 

6320 ± 10 
6340 

6330 ± 20 



6326±f 



6194 ~ 6292 6284 ~ 6357 65 < A 1S < 90 
6253TJ 3 632811 69 < A 1S < 80 



62581? 



ii 



63331? 



ii 



* The experimental mass of this singlet state is presented in [2]. 
t Averaging over the five values in Table I. 

* We treat Eichten and Quigg's results [1] in the same manner. 
T Our best estimation for the center-of-gravity triplet state. 
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TABLE III. Parameters used for Igi-Ono potential. 



A M (GeV) a {GeV 2 ) g {GeV) d {GeV 2 ) m c (GeV) m b (GeV) 



o.it 


0.1733 


.3076 


0.4344 


1.134 


4.563 


0.2 


0.1587 


0.3436 


0.2550 


1.322 


4.731 


0.3 


0.1443 


0.3280 


0.0495 


1.471 


4.868 


0.4 


0.1387 


2.903 


0.582 


1.515 


4.910 


0.5 


0.1391 


2.955 


1.476 


1.514 


4.911 


0.1^ 


0.1762 


0.2753 


0.4720 


1.120 


4.551 


0.2 


0.1734 


0.3479 


0.5362 


1.267 


4.684 


0.3 


0.1615 


0.4482 


0.6020 


1.416 


4.815 


0.4 


0.1389 


0.6219 


0.5632 


1.604 


4.986 


0.5 


0.1137 


1.0029 


0.7368 


1.748 


4.118 



' Type I potential, 
t Type II potential. 
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TABLE IV. The cb mass spectra predicted for various 
using Igi-Ono (type I and II) potential (all in MeV). 



States Ref.[6,24] A^ 


= 100 


200 


300 


400 


500 


b = 


20 












IS 


6327 


6355 


6342 


6336 


6347 


6349 


2S 


6906 


6941 


6928 


6907 


6911 


6923 


3S 


7246 


7290 


7266 


7270 


7272 


7274 


IP 


6754 


6781 


6768 


6759 


6763 


6751 


2P 


7154 


7170 


7155 


7151 


7153 


7143 


ID 


7028 


7055 


7041 


7030 


7031 


7019 


2D 


7367 


7360 


7351 


7353 


7354 


7342 


b = 


5 












IS 


6327 


6357 


6347 


6342 


6326 


6316 


2S 


6906 


6940 


6921 


6936* 


6937 


6930 


3S 


7246 


7284 


7300 


7262 


7272* 7238 


IP 


6754 


6782 


6766 


6763 


6749 


6746 


2P 


7154 


7168 


7161 


7160 


7139 


7136 


ID 


7028 


7055 


7038 


7038* 


7026 


7022 


2D 


7367 


7361 


7346 


7340 


7335 


7347* 



* Carried out to the first order. 

* Carried out up to the second order correction. 
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TABLE V. The cb mass spectra and A n 5 for various Ajjg 



calculated by using Igi-Ono potential with a s = 0.250 (all in MeV). 
b= 20 5 

States Ajfs = 100 200 300 400 500 100 200 300 400 500 
IS 6327* 6329 6316 6310 6321 6323 6331 6321 6316 6300 6290 
1% 6276 6263 6256 6268 6273 6277 6271 6265 6250 6239 

l 3 Si 6347 6334 6328 6339 6339 6349 6338 6332 6317 6307 

Ais 71 71 72 72 67 72 67 67 67 68 

2S 6906* 6915 6902 6881 6885 6897 6914 6895 6910 6911 6904 
2 1 S 6888 6876 6856 6861 6873 6888 6869 6883 6883 6878 

2 3 S 1 6924 6911 6889 6893 6905 6923 6904 6920 6920 6912 

A 2S 36 35 33 32 32 35 35 37 37 34 



*Here we cite Ref. [6]. 
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TABLE VI. The cb mass spectra and A nS using 
Igi-Ono (type III) and Chen-Kuang potentials. 



A-MS — 


300 


300 


390 


100 


180 


375 


b = 


16.3 


19 


19 








m c = 


1.506* 






1.478 f 






m b = 


4.897* 






4.876* 






M(cb) 














IS 


6309* 


6337 


6298 


6323 


6323 


6323 


l 3 ^ 


6326 


6354 


6318 


6340 


6340 


6340 


1% 


6258 


6287 


6238 


6272 


6272 


6272 


Ais 


67 


67 


81 


68 


68 


68 


2S 


6880 


6898 


6878 


6879 


6879 


6879 


2 3 S 1 


6889 


6906 


6886 


6888 


6888 


6888 


2^0 


6855 


6874 


6853 


6853 


6853 


6853 


A 2S 


33 


33 


33 


35 


35 


35 


35 


7247 


7262 


7255 


7257 


7257 


7257 


45 


7553 


7569 


7564 


7570 


7570 


7570 


IP 


6725 


6749 


6738 


6723 


6723 


6723 


2P 


7124 


7142 


7136 


7126 


7126 


7126 


ID 


6997 


7018 


7014 


6992 


6992 


6992 


2D 


7328 


7345 


7342 


7332 


7332 


7332 



*.Mass (in GeV) fitted for Igi-Ono (type III) potential, 
t Mass (in GeV) fitted for Chen-Kuang potential. 
* Here we cite Ref. [6]. 
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TABLE VII. The fine and hyperfine splittings in our work compared with 
that in other works.. 



States Quadratic* Linear* Fulcher* Lattice* C-Q t I-O* 



b = 










16.3 










100- 


375 300 


Fine splittings 












M(2S)-M(1S) 


558* 


533 


579 


672 ± 120 556 


571 


M (35) - M(1S) 


931 


899 




934 


938 


M(2P) - M(IP) 


382 


376 




403 


399 


M(2D) - M(ID) 


324 


321 




340 


331 


b = 




20 f 


5t 


19 1 " 19 f 


16.3+ 






400 


100 


390 300 


300 


M(2S)-M(1S) 


558 


564 


583 


580 561 


571 


M(3S)-M(1S) 


931 


925 


927 


957 927 


938 


Hyperfine splittings 












Ais 


76 


62 


55 


41 ± 20 68 


67 


A 2S 


42 


33 


32 


30 ± 8 35 


33 



* Here we cite Ref . [6] . 
t We used SLNET. 



* Splitting masses (in MeV).. 
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